function var=dLdDelta(t,y,rho,kW,RefWealth,sigma,piW,n)
var=y*(-piW-rho+kW*(t-RefWealth)^(-sigma)*y)/((-piW-n)*t);